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Abstract 



In this paper, a compensated compactness framework is established for sonic-subsonic ap- 
proximate solutions to the ri-dimensional(n > 2) Euler equations for steady irrotational flow 
^ ■ that may contain stagnation points. This compactness framework holds provided that the ap- 

\ proximate solutions are uniformly bounded and satisfy 7?;~^(ri) compactness conditions. As 

illustration, we show the existence of sonic-subsonic weak solution to n-dimcnsional(n > 2) 
I Euler equations for steady irrotational flow past obstacles or through an infinitely long nozzle. 

This is the first result concerning the sonic-subsonic limit for n-dimension(rt > 3). 

■ Key wordsMulti-dimension, Sonic-subsonic flow. Steady irrotational flow. Compensated com- 

\ pactness 
^ ■ 2000 MR Subject Classification 35L65 

1 Introduction 

p\j ' The n-dimensional(n > 2) Euler equations for the steady irrotational flow reads 

{curl u = 0, 
div(H = 0, (1.1) 
div(/3u ®u + pi) = 0, 



where u = (ui, • • • are the flow velocities, and p and p represent the density and pressure- 
density function, respectively, and (curl u)ij = dx^Ui — dx^Uj, i, j = 1, ■ ■ ■ ,n, is a n x n matrix 
and J is a n X n unit matrix. Usually, we require p'{p) > for p > 0. 

The famous Bernoulli's law can be easily derived from (jl.ip i and (jl.iP '^: 



Hp) + 2<1 = const, (1.2) 



where = = Yl^=i '^1 the flow speed and h{p) = ^^-j^ds is the enthalpy. 

In this paper, we are interested in the po^ 
(11.21) is converted to the normalized formula 



In this paper, we are interested in the polytropic gas, that is p = p{p) = for 7 > 1. Then 



7 — 1 n. ^ 

The local sound speed is defined by 



p = p{q) = {1 - l-±q^)^ . (1.3) 



„2_ 7-l„2 



p'ip) = 1 - ^q'- (1-4) 
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At the sonic point q = c, ()1.4p implies = The critical speed qcr is defined as 




and the Bernoulli's law is rewritten as 



Thus the flow is subsonic when q < Qcr, sonic when q = qcr, and supersonic when q > Qcr- 

For the isothermal flow, p = c^p, where c > is the constant sound speed, Bernoulli's law 
takes: 



for some constant po > 0. 

It is well known that the steady irrotational Euler equations (jl.ip is of mixed type of partial 
differential equations, which is elliptic if g < qcr, parabolic if q = qcr, hyperbolic if q > Qcr- Since 
the equations of uniform subsonic flow possess ellipticity, its solutions have extra-smoothness 
to those related to transonic flow or supersonic flow. There are a large of literatures on the 
smooth uniform subsonic solutions, for instance, see [5] . |29j . [32] . |47] . [38] for two dimensional 
flow and [M] , [25] , [26] , [30] , [H] , [51] , [52] , [53j for three dimensional flow. Among them, Frankl and 
Keldysh [32] obtained the first result about the subsonic flow past a two dimensional flnite body 
(or airfoil). Shiffman in [47] . |48j proved there exists a unique subsonic potential flow around 
a given proflle with flnite energy provided that the infinite free stream fiow speed q^Q is less 
than some critical speed, which was improved by Bers [5]. Finn and Gilbarg [29j proved the 
uniqueness of the two dimensional potential subsonic flow past a bounded obstacle with given 
circulation and velocity at inflnity. The flrst result for three dimensional subsonic flow past an 
obstacle was given by Finn and Gilbarg [30] in which they studied the existence, uniqueness 
and the asymptotic behavior with some restriction on Mach number. Dong [23] extended the 
results of Finn and Gilbarg [30] to maximum Mach number M < 1 and to arbitrary dimensions. 
Recently, Du, Xin and Yan [26] obtained the smooth uniform subsonic for n-dimensional(n > 2) 
flow in an inflnitely long nozzle. For other related results, we refer to [2l [3l HI [6l [28 t [3 H [33 t [33 1 [35] 
and references therein. 

However, few result is known until now for the sonic-subsonic flow and transonic flow, because 
the uniform ellipticity is lost and shocks may present. That is, smooth solutions may not 
exist, and weak solution is necessarily considered. Morawetz \A0\ [3T] firstly introduced the 
compensated compactness method to study steady fiow of irrotational Euler equations. Indeed, 
Morawetz established a compactness framework under assumption that the stagnation points 
and cavitation points are excluded. Morawetz's result was improved by Chen, Slemrod and 
Wang [15j in which the approximate solutions are constructed by a viscous perturbation. 

On the other hand, the first compactness framework on sonic-subsonic irrotational flow in 
two dimension was recently due to Chen, Dafermos, Slemrod and Wang [14J and Xie-Xin [51] 
independently by combining the mass conservation, momentum, and irrotational equations. 
The key point of [13] and [51] is based on the fact that the two dimensional steady flow can be 
regarded as the one dimensional system of conservation laws, that is, x is regarded as time t, so 
that the div-curl lemma can be applied to the two momentum equations. In fact, the authors 
|14j flrst applied the momentum equations to reduce the support of the corresponding Young 
measure to two points, then the irrotational equation and the mass equation are used to deduce 
the Young measure to a Dirac measure. Xie and Xin [52j investigated the sonic-subsonic limit 




(1.5) 
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for the three-dimensional axis-symmetric flow(it is similar to the two dimensional case) through 
an axis-symmetric nozzle. 

However, the compactness framework established in [T3] is no longer effective for n-d (n > 3) 
steady irrotational Euler equations, which can not be reduced to one dimensional system of 
conservation laws, and the famous div-curl lemma is no longer valid for the momentum equations. 
In this paper, we find that it is enough, by only using the mass conservation equation and 
irrotational equations, to reduce the Young measure to a Dirac measure for arbitrary dimension. 
Thus we establish a compactness framework of approximate solutions for steady irrotational 
flow in n-dimension(n > 2). It is worthy to point out that the famous Bernoulli's law plays a 
key role in our proof. As application, we show the sonic-subsonic limit for steady irrotational 
flows past obstacles or through an infinitely long nozzle in n-dimension(n > 2). 

The rest of this paper is organized as follows. In section 2, we establish the compact- 
ness framework of sonic-subsonic approximate solutions for the system of steady irrotational 
equations in n-dimension(n > 2). In section 3, we give two applications of the compactness 
framework to show the existence of sonic-subsonic flow over obstacles or through an infinitely 
long nozzle. Finally, in section 4, we extend the compensated compactness framework to the 
realm of self-similar solutions to the Euler equations for unsteady irrotational flow. 

2 Compensated Compactness Framework for Steady Irrotaitional 
Flow in n-Dimension 

Let a sequence of function u^{x) = {uf, ■ ■ ■ ,u^){x), defined on open subset 0, C M", satisfy the 
following Conditions: 

(A.l) q^{x) = \u^{x)\ < Qcr a-e. in Q; 

(A. 2) curl u^, and div(p(g^)M^) are confined in a compact set in H^^{n). 

Based on the above conditions, the famous div-curl lemma and the Young measure repre- 
sentation theorem for a uniformly bounded sequence of functions imply: 

n 

< p{Q)q^, v{u) >= < ^{u) X Ui, v{u) >, (2.1) 

i=l 

where v = Vx{u) is the associated Young measure (a probability measure) for the sequence 
u'^{x) = (uf,--- ,ti^)(x). Now, the main effort is to establish a compensated compactness 
framework, namely, to prove that is a Diac measure by using the identity ()2.ip . This in turn 
implies the compactness of the sequence u^{x) = (uf , • • • ,«^)(x) in L\^^{Q). 

Theorem 2.1 (Compensated compactness framework) Let a sequence of functions u^{x) = 
(nf,--- ,ti^)(2;) satisfy conditions {A.l) and (^.2). Then the associated Young measure u is 
a Dirac mass and the sequence li^(x) is compact in Lj^^{Q,); that is, there is a subsequence (still 
labeled) u^{x) u{x) = {ui, • • • ,Un){x) a.e. as e — )■ and satisfying q{x) = \u{x)\ < qcr, a.e. 
X eQ. 

Proof. Let 

n 

/(^.«,n(^)) = -nf))(p(,«)n« -p(,(2))^f)), (2.2) 

1=1 

where n*-*^ = {ui \ ■ ■ ■ ,Un^) and g^*^ = |n*-*^| for z = 1,2 be two independent vector variables. 
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After some basic calculations on I{u^^\u^'^^), we have 

n n 
i=l i=l 

Then the Cauchy inequality implies 

= (2))2rfM(,~). (2.4) 



where q lies between (7^^^ and q^"^^ due to the mean-value theorem. The famous Bernoulli's law 
(|1.3|) gives that for 7 > 1, 

which immediately implies 



_ ,1 _ 1^,2)^-1(1 _ l±i,2) ^ ,1 _ 1^,2)^-'(1 - (2.5) 



For 7 = 1, the Bernoulli's law is 

p{q) = po exp 

which gives 



2cr 



2^7' 



_.p„(l__)exp(^-^j. (2.6) 

Thus , for sonic-subsonic flows, namely, g^^^g^^^ < gcr, (|2.4p — (|2.6p imply 

/(^z«,n(2)) = (^(1) _g(2))2^(g~) >o (2.7) 

and 

(g(i) _ ^{2))2^(^~) ^ 0^ if only if g(^) = (2.8) 



From the identity ()2.ip . noticing that the Young measure is a probability measure, we have 

</(n«,u(2)),i.(n«)®z.(n(2)) >= 0, (2.9) 

which together with ()2.7p and ()2.8p implies f;^"^^ = q^'^\ where v{u^^^) iy{u^'^^) is understood as 
a product measure of u{u^^^) and z^(ti(2)). "With the property q^^^ = q^"^^ at hand, we have from 

(1221) 

= </(nW,'u(2)),i.(nW)^i/(n(2)) > 

n 

= < PiQ^'^) - ^^(^^'^) ^ H^^^^^) >, (2.10) 

1=1 

which immediately implies , the Young measure is a Dirac measure. This completes 

Theorem 12. li 



5 



Remark 2.2 Theorem \2.1\ is valid for any n > 2. Namely, a compactness framework in The- 
orem \2.1\ is established for sonic-subsonic limit for steady irrotational flow in arbitrary dimen- 
sion. From the Bernoulli's law ()1.2p . it is straightforward to extend Theorem \2.1\ to the general 
pressure- density function p satisfying p'{p) > for p > 0. 

Theorem \2.1\ is valid for any n > 2. Namely, a compactness framework in Theorem \2.1\ is 
established for sonic-subsonic limit for steady irrotational flow in arbitrary dimension. From the 
Bernoulli's law ()1.2p . it is straightforward to extend Theorem \2.1\ to the general pressure- density 
function p satisfying p'{p) > for p > 0. Futhermore, The above idea can be applied for the 
sonic-subsonic limit of the self- similarity solutions, that is, {p,u){x,t) = {p,u){j), x € i?", for 
unsteady compressible Euler system. The idea is also available for Euler-Poisson systems. 



We now consider a sequence of approximate solutions to the Euler equations (jl.ip i . 
(jl.ip o and the Bernoulli's law (|1.3p or ()1.5p . That is, besides Conditions (A.l) and {A.2), the 
approximate solutions further satisfy 

{curl u'^ = oiie), , , 

^ " (2.11) 
div(/)(g^)n^) = 02(e), 

where oi(e), 02(e) — > in the sense of distributions as e — > 0. Then, as a corollary of Theorem 
12.11 we have 

Theorem 2.3 (Convergence of approximate solutions) Let u^{x) = (uf,-- - ,u^){x) be a se- 
quence of approximate solutions satisfying ()2.1ip and the Bernoulli 's law ()1.3p or ()1.5p . Then, 
there is a subsequence (still labeled) u^(x) that converges a.e. as e to a weak solution 
u{x) = («!,■■■ ,Un){x) to the Euler equations of p.ip i . (|l.ip 9 and the Bernoulli's law (jl.3p or 
(jl.Sp satisfying q{x) = \u{x)\ < qcr, a-e. x G Q. 



Remark 2.4 For any functions Q{u) = (Qi(u), • • • , Qn{u)) satisfying 

div{Q{u^)) = o(e), (2.12) 

where o(e) — >■ in the sense of distributions as e 0, from the strong convergence of , 
Q{u) = holds in the sense of distributions. So if we have 

div{p{q^)u^ (8) + p'^I) = o(e) 0, in the sense of distributions, (2-13) 

the weak solution in Theorem 12.31 also satisfies the momentum equations (jl.ip q in the sense of 
distributions. 



There are various ways to construct approximate solutions by either numerical methods or 
analytical methods such as vanishing viscosity methods. In the next section, we will show two 
examples to apply the compactness framework built in Theorem 12.11 

3 Sonic Limit of Irrotational Subsonic Flows in n-Dimension 

In this section, we wish to apply the compactness framework established in Theorem 12.11 to 
obtain the sonic limit of n-dimensional(n > 2) steady irrotational subsonic flows. 

Firstly, we give an example of subsonic-sonic limit past obstacles. Let the obstacle T be one 
or several closed n — l(n > 2) dimensional hypersurfaces. We shall always assume F G (7^'"^°. 
Denote by I'(F) the domain exterior to F, see Fig 13. H 
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Figure 3.1: General high dimensional case. 



Problem P(noo) : Let n > 2. Find functions u = {ui, ■ ■ ■ satisfy (jl.ip i . ()l.ip 9 with the 
Bernoulli's law (jl.3p or ()1.5p . and the slip boundary condition 

(pu) ■ n = on r, (3.1) 

where n denotes the unit inward normal of domain T>{T), and the limit 

Uoo = lim u{x), (3.2) 

exists and is finite. 

The main result of [23] is described as follows: 

Theorem 3.1 (Uniform Subsonic Flows Past An Obstacle for n-D Case |24j ) Let 

QoQ '■= \uoo\- There exists a positive number q < qcr, so that P(uoo) has a uniform subsonic 
solutions if < Qoo < q- Furthermore, let qm{qoo) = sup |n(x)|, then the function qmiqoo) £ 

C[0, q) and qm{qoo) Qcr as q^o q- 

We note that the above Theorem does not apply to the critical flows, that is, those flows for 
which qmiqoo) = qcr and which hence must be sonic at some point. Here, by Theorem 12. H we 
establish a more general result. 

Theorem 3.2 (Sonic Limit Past An Obstacle) Let — )• u^o as e ^ be a sequence of 
speeds at oo with q^ < q = \uoo\, and = (nf,--- ,u^) be the corresponding solutions to 
Problem P(n^). Then, as Uoo, the solution sequence u^{x) possess a subsequence (still 

denoted by) u^{x) that converge a.e. in V(T) to a vector function u{x) = (ui, • • • , Un){x) which is 
a weak solution of Problem P(noo) with goo = q- Furthermore the limit velocity u = {ui, ■ ■ ■ ,Un) 
also satisfies (jl.ip ci in the sense of distributions and the boundary conditions ()3.ip as the normal 
trace of the divergence-measure field (ui, ■ ■ ■ ,Un) on the boundary (see fllf)- 

Proof. The strong solutions satisfy (jl.ip . and the Bernoulli's law and are uniform subsonic 
solutions of Problem P(u^). Hence Theorem 12.11 immediately implies that the Young measure is 
a Dirac mass and the convergence is strong a.e. in T>{T). The boundary conditions are satisfied 
for u in the sense of Chen-Frid [llj. On the other hand, Since (jl.ip ci holds for the sequence of 
subsonic solutions u^{x), it is straightforward to see that u also satisfies (jl.ip -? in the sense of 
distributions. This completes the proof of Theorem 13.21 
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Now we give another example of subsonic-sonic limit through an infinite long nozzle. As in 
[26| . denote the multi-dimensional nozzle domain by 0, which satisfies the following regularity 
assumption: there exists an invertible C^'" map T : Q ^ C : x ^ y satisfying 

'T{dn) = dC, 

For any /c G M, T{n fl {x„ = k}) = B{0, 1) x {y^ = k}, 
< \\T\\c2,c,\\T-'^\\c2,c. < K < oo, (3.3) 
The nozzle approaches to a cylinder in the far fields, i.e, 
^1 n {xn = k} ^ S± as k ^ iboo, respectively, 

where K is a uniform constant, C = -6(0,1) x (— oo,+oo) is a unit cylinder in R", i?(0, 1) is 
unit ball in M""^ centered at the origin, S± are n — 1 dimensional simply connected C"^'", x„ is 
the longitudinal coordinate and x' = {xi, ■ ■ ■ ,x„_i) G M""-*^, see Fig l3.2t 




Figure 3.2: n-dimensional nozzle 

Problem P(?7io) : Let n > 2. Find functions u = (ui, • • • ,Un) satisfy (jl.ip i . (jl.ip 9 with the 
Bernoulli's law ()1.3p or (jl.Sp . and the slip boundary condition 

{pu) • n = on dn, (3.4) 

where n denotes the unit outward normal of domain 0,; and the mass flux condition 

/ p{\u\^)u ■ tdS = mo > (3.5) 
Jso 

The main result of [26j is stated as follows: 

Theorem 3.3 (Uniform Subsonic Flows in n-D Nozzle[26j) There is a critical mass flux 
Mc > 0, which depends only on Q, such that ifO < mo < Mc, then P(mo) has a unique uniformly 
subsonic flow through the nozzle, i.e, qmif^o) := sup|n(x)| < Qcr- The velocity (ui,--- ,Un) is 

Holder continuous. Moreover, Qmin^o) ~^ Qcr cls uiq Ale. 

Similar to Theorem 13.21 we have 

Theorem 3.4 (Sonic Limit Through A Nozzle) Let < tuq < Mc be a sequence of mass 
flux, and let u^{x) be the corresponding solution to P(mQ). Then as ttIq Mc, the solution 
sequence u'^{x) possess a subsequence (still denoted by) u^{x) that converge strongly a.e. in Q 
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to a vector function u{x) = (ui, • • • ,Un){x) which is a weak solution ofP{Mc) with Bernoulli's 
law. Furthermore the limit velocity u = (ui, ■ ■ ■ ,Un) also satisfies in the sense of distri- 

butions and the boundary conditions ()3.4p as the normal trace of the divergence-measure field 
(ui,--- ,Un) on the boundary (see flT] / ). 

Remark 3.5 In this section we only give two examples as applications of compactness frame- 
work established in Theorem \2.1\ Certainly it can be used to other cases as long as conditions 
{A.l) and {A.2) are satisfied. 
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